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Abstract 

This article presents the Serial and Unserial Methods (SUM). The algorithms are 
strongly related to the first part of a classical reference in combinatorics, the Com¬ 
binatorial algorithms for computers and calculators, from Albert Nijenhuis and Her¬ 
bert Wilf. The Serial Method proposal is to obtain the output of a specihc kind of 
combinatorial family from its position on the list of all combinatorial possibilities. 
The Unserial Method is the inverted step of Serial Method, getting the serial num¬ 
ber from the combinatorial family given as input. The serial number is the position 
of the combinatorial family on the list. 
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1 Introduction 


On the classical reference of combinatorial algorithms [1], A. Nijenhnis and 
H. Wilf presents the Next algorithms, and the proposal of this methods is 
to obtain the next ontpnt of a combinatorial family from the actnal one. 
Getting an example from Nijenhuis-Wilf’s reference, the Next Permutation [1] 
algorithm generates de whole list of permntations withont considering the list 
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position of each permutation contained in the list. This is a fast way to obtain 
all the permutations or a set of permutations from a specihc one. This work 
is a set of algorithms which are concentrated to get a specihc combinatorial 
family on the list. The Serial and Unserial combinatorial family in this article 
is the permutation, composition, partition of an n-set, /c-subsets and subsets. 
The combinatorial family exposed in this article is divided on the sections 
below. 


2 Serial Permutation Method (SPM) 

2.1 Basic Concepts 


The permutation algorithms are classihed in two main groups: one that creates 
a set of permutations from the identity permutation and other that produces a 
set of permutations by means of simple changes between the vector elements, 
creating a new permutation from the actual one. The algorithm Next Permu¬ 
tation for N letters [1] belongs to the second group, creating a complete set 
of permutations through successive algorithm invocations, getting the next 
permutation until reaches the last vector of the list. The algorithm Next Per¬ 
mutation is creates the next permutation vector using only local information. 
The next vector is different from anyone else that was dehned before. 

There is a specihc problem on this creating process: 

Is it possible to get directly a permutation vector located on a specific position 
on the list of vectors, excluding the alternative of the Next Permutation 
algorithm’s successive execution until gets the desired vector? 

The Serial Permutation Method (SPM) answers this question, being able to 
process the same list of permutation that the Next Permutation algorithm 
does, with a difference: the SPM only needs the serial number and the permu¬ 
tation vector’s size to process the desired output. It is also possible to invert 
this process from the SPM, getting the serial number through the permutation 
vector. 


2.2 SPM’s Construction 


The SPM was developed from the observation of the Next Permutation’s data 
outputs. This algorithm has an auxiliary variable called offset vector (also 
called inversion vector) [1], which is dehned below. 
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Consider p the permutation vector with n elements and let d be the offset 
vector with n — 1 elements, the d vector is defined by [1]: 

di = \j ■■ j < i, Pj > Pi+i\ i = l,2,...,n (1) 


That is explained: 

How many numbers are there biggest than pj+i between the beginning of the 
permutation vector and the i index of the d vector? 

Being n\ the complete list of the created vectors, we have in Table 1 an example 
of output from the p permutation vector related to its serial, for n = 5 and 
the serial = 32. 

Table 1 

Input and output processing that are common in the Next Permutation algorithm. 


Input 

Output 

n = 5 

p =(3,5,1,2,4) 

Serial = 32 

d= (0,2, 2,1) 


The SPM is subdivided in two steps, and the data input of the second step 
corresponds to the data output of the first step: 

(1) Given a serial number and the size of the permutation vector, determine 
the offset vector; 

(2) Given the offset vector, determine the permutation vector. 

2.3 Algorithm to Attain the Offset Vector 


To solve the first SPM step it is necessary an algorithm to get the offset vector, 
which is called Serial Offset Algorithm (SOA). The algorithm does a calculus 
that reflects the pattern of the offset vectors’ creation throughout the whole 
list of permutation, using only local information. In Table 2, we can visualize 
the pattern of the offset vector in the whole list of permutation, where the 
detached column will serve as a guide to the determination that the SOA 
decodes. 

Let d be the offset vector with n elements. Each d vector’s column showed in 
Table 2 has the following properties: 

• Let k be the value of the current column {k = 1, 2,..., n), the elements of the 
dk column (0 < < k) repeat k\ times until the element of the column gets 

the k value. This can be visualized on the first part of the detached column 
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Table 2 

Full list of permutations, with n = 4 


Serial 

P 

d 

1 

(1,2,3,4) 

(0,0,0) 

2 

(2,1,3,4) 

(1,0,0) 

3 

(3,1,2,4) 

(1,1,0) 

4 

(1,3,2,4) 

(0,1,0) 

5 

(2,3,1,4) 

(0,2,0) 

6 

(3,2,1,4) 

(1,2,0) 

7 

(4,2,1,3) 

(1,2,1) 

8 

(2,4,1,3) 

(0,2,1) 

9 

(1,4,2,3) 

(0,1,1) 

10 

(4,1,2,3) 

(1,1,1) 

11 

(2,1,4,3) 

(1,0,1) 

12 

(1,2,4,3) 

(0,0,1) 

13 

(1,3,4,2) 

(0,0,2) 

14 

(3,1,4,2) 

(1,0,2) 

15 

(4,1,3,2) 

(1,1,2) 

16 

(1,4,3,2) 

(0,1,2) 

17 

(3,4,1,2) 

(0,2,2) 

18 

(4,3,1,2) 

(1,2,2) 

19 

(4,3,2,1) 

(1,2,3) 

20 

(3,4,2,1) 

(0,2,3) 

21 

(2,4,3,1) 

(0,1,3) 

22 

(4,2,3,1) 

(1,1,3) 

23 

(3,2,4,1) 

(1,0,3) 

24 

(2,3,4,1) 

(0,0,3) 


in Table 2, illustrated in blue, which the number of detached elements in 
blue or red is dehned by {k + 1)! 

• When dk = k in [k + 1)\ times, the following elements from the list are put 
on the inverted form, like showed in red in Table 2. The direct and inverted 
list intercalate themselves until complete all the positions of the column. 

Exists throughout each offset vector column an intercalation between the two 


4 






























kinds of lists (direct and inverted), we can consider this fact as an element of 
parity in the list, considering as an even the list in its direct form and as odd 
the list in its inverted form. To determinate the relation between the serial 
number and the parity of the list, follows: 


/ 


s — 1 
.(fc + 1)!. 


mod 2 


( 2 ) 


Where: 

• /: Determines if the list is direct (/ = 0) or inverted (/ = !); 

• s: Serial number; 

• k: Index of the offset vector; 

• |_xj: Floor function. It returns the biggest integer value smaller than x] 

• mod: An operation that returns a division’s rest. 

If the list had been direct, we attribute to the elements of the offset vector: 
[{s — 1) mod {k + 1)\] 


If the list had been inverted, we attribute the complement which would be the 
direct list: 


dk 


k- 


[(s — 1) mod {k + 1)!] 
k\ 


( 4 ) 


Each attribution is made through a loop that go through the offset vector. 
The equations 3 and 4 can be joined in a unique form, but is a inefficient way 
to compute the equation: 


dk 


f.k+ i-iy. 


[(s — 1) mod {k + 1)!] 
k\ 


( 5 ) 


Where / is defined on equation 2. 

Serial Offset Algorithm (SOA) 

Routine Specifications: 

• n: Size of the offset vector; 

• i: Index of the offset vector; 

• d: Offset vector, alternating its indices from 0..n — 1; 

• s: Permutation’s serial. 

Routine: 
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for i 1 to n do 

if [(s — !)/(* + 1)!J mod 2 = 1 then 
di-i ^ i - [((s - 1) mod {i + 1)!) / i!)J 
else 

(ij_i ^ [((s — 1) mod {i + 1)!) / i!)J 
end for 
return d. 


2-4 Algorithm to attain the permutation vector from the offset vector 


After the SOA is computed, the SPM is concluded with the Permutation Al¬ 
gorithm by Offset (PAO). The PAO does the SPM’s second step, returning the 
desired output. A relevant topic for the PAO construction is to hud the de¬ 
coding process of the permutation vector, with only the offset vector being the 
input. We know that the offset vector maps the elements of the permutation 
vector. In the decoding, we have: 

Let p be the permutation vector and d the offset vector: 


P = (Pl,P2, •••,Pn) d = {di,d2, ...,dn-l) 


( 6 ) 


And d already has its values determining by the SOA. This is how the decoding 
process is made: we know that di has its values discretely included between 0 
and 1. Obviously, 0 and 1 are the unique possible elements for di. From this 
information, we can conclude: 


di = 


0, if Pi < p2 
1, if P2 < Pi 


We do not consider the possibility of equality between the elements of the 
permutation vector, because we know that there is no repeated elements on 
the vector, and if we organize it, the difference between them will be only one 
unit. Like di, the element d 2 has its values included between 0 and 2. In this 
and in other cases, we analyze the current inequality from the inequality that 
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was created previously. For ^ 2 , we have: 


0 , d2 — 


0, if Pi < P2 < P3 
< 1, if Pi < P3 < P2 
2, if P3 < Pi < P2 


di = < 



0, if P2 < Pi < P3 
< 1, if P2 < P3 < Pi 
2, if P3 < P2 < Pi 


We can conclude that as the value of ^2 increases, p 3 ’’slides” on the left 
through the inequality. In general, we have: 

4 e {0, l,...,/c} {pi < pj < ... < pfc < ... < pj (7) 

With pfc in the inequality with dk positions, counted from right to left, because 
the order of the elements is ascendent. Finished the offset vector’s raster, we 
will get a set of inequalities that informs the order of the permutation elements. 
We have, for example: 

Given a d = (0, 2, 2) offset vector, the inequality for a p permutation vector 

is: 


P3 < P4 < Pi < P2 


( 8 ) 


Being the last inequality at 8 the hnal disposition between the elements of the 
vector. After the mapping of the permutation’s elements was done, we can say 
that each element of the vector ordered in 8 corresponds to the elements of the 
identity permutation. This fact classihes the SPM in the hrst group described 
in the basic concepts of this section. On the 8 inequality, we have: 

(P3 = 1) < (P4 = 2) < (pi = 3) < (P2 = 4) (9) 

The next step to get the hnal output is to arrange each one of the elements 
ordered on its own positions. As the p vector is ordered like (pi,P 2 , ■■■,Pn), we 
have: 

P = (Pi,P2,P3,P4) ^ p=(3,4,1,2) (10) 

So we can get the SPM’s hnal output. In implementation terms, to compute 
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the input built from the permutation vector (from right to left) the corre¬ 
sponding indices of the vector were used in relation to its complement. We 
will notice at the implementation that in the insertion moment of the cur¬ 
rent element, if it replaces another, the elements of inequality located on the 
left will ’’slide” to the left side, allocating space for the current element. The 
’’slides to the left” operation is implemented on the subroutine push. 

Permutation Algorithm by Offset (PAO) 

Routine Specihcations: 

• n\ Size of the permutation vector; 

• i,j\ Indices of the algorithm’s vectors; 

• p: Permutation vector, alternating its elements on 0..?7, — 1; 

• r: Vector which will keep the element’s position before being ordered. 


Routine: 



1 


if Tn-l-di.i = 0 


1 to n — 1 do 
= 0 then 



else 


<— push 


end for 


for i <— 0 to n — 1 do 



end for 
return p. 


Subroutine Push 

for j <— (n — 1) — i to j < (n — 1) — (ij_i do 


^ 0+1 

end for 
return i + 1. 


3 The Inverted Process of SPM (Unserial Method) 


We can describe now the inverted process of the SPM’s. We have the permu¬ 
tation vector as the input, and the desired output is the correspondent serial 
number. The ingenuous process to get the serial number for the permutation 


vector is the raster of the permutation list, comparing the vectors one by one, 
until gets the equivalent vector computed on the input, being the returned 
value the loop’s index that does this raster. However, we can hnd the serial 
value inverting the SPM’s steps: 

(1) Given the permutation vector, hnd the offset vector; 

(2) Given the offset vector, hnd the serial number. 

Like the SPM, the second process depends on the hrst, with the hrst step data 
output corresponding to the second step data input. 


3.1 Algorithm to attain the offset vector from the permutation vector 


We will call this algorithm as the Offset Algorithm by Permutation (GAP). 
It uses on details the dehnition of the ohset vector [1] (see equation 1). Two 
nested loops add the value of each element of the ohset vector. 

Offset Algorithm by Permutation (OAP) 

Routine Specihcations: 

• n: Size of the permutation vector; 

• i,j: Indices of algorithms’s vectors; 

• p\ Permutation vector alternating its elements on 0..n — 1; 

• d\ Ohset vector alternating its elements on 0..n — 2. 

Routine: 

for i 0 to n — 2 do 
for j ^ 0 to i do 

if pj > Pi+i then 
di ^— dj + 1 

end for 
end for 
return d. 


3.2 Algorithm to attain the serial number from the offset vector 

Like the OAP, the Serial Algorithm by Offset (SAO) also uses dehnitions 
already showed at this section. We know that the value attributed for the 
elements of the ohset vector, dehned on 3 and 4 can be ordered depending on 
the serial. Given the variables: 
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• d = ((ii, (^ 2 ,(ifc, dn)'- Offset vector; 

• dk- Offset vector’s element; 

• s: Permutation’s serial number; 

• k: Index of the offset vector; 

• q: Quotient of the division between s — 1 and {k + 1)!; 

• |_a;J: It returns the biggest integer value smaller than x. 

For direct list, we have: 

s - 1 = dk-kl + [q\.{k + 1)\ (11) 


For inverted lists, we have: 

s - 1 = {k - dk).k\ + [q\.{k + 1)! 


( 12 ) 


An important question for the algorithm implementation is how to hnd the 
variable’s q value. To solve this question we have to consider that dn belongs 
to a column that has only one list on a direct disposition (see Table 2, last d 
vector column). With this information, we conclude that the quotient for this 
column corresponds to a value between [0,1 [, making: 

L«J.(fc+l)! (13) 


Corresponds to zero. With this, 11 and 12 is equivalent to: 

s — 1 = dk-k\ (14) 

For direct lists, and 

s — 1 = {k — dk)-k\ (15) 


For inverted lists. 

The strategy of implementation to get the serial is done by incremental mode, 
working the current information being based on past information, doing a 
raster on the offset vector from right to left. This is particularly useful on 
this case, considering the next value that will converge to the final serial that 
belongs to a group related with the previous serial elements. Each quotient 
q during the raster will corresponds to zero, because we are considering the 
division in relation to the gotten serial. Another relevant question is how to 
classify if the previous element of the offset vector it is contained in a direct 
or inverted list. Looking at Table 2, we can make easily an equivalence table, 
showed in Table 3. 
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Table 3 


Next element status list o 

’ the offset vector 

Previous Value - Parity 

Next Value - Parity 

Element 

List 

List 

Even 

Direct 

Direct 

Even 

Inverted 

=> Inverted 

Odd 

Direct 

Inverted 

Odd 

Inverted 

=> Direct 


Exemplifying the attainment serial process from the offset vector, we have: 

Let d = (1, 0, 3) be the offset vector and s the permutation serial. As the raster 
is done from right to left, we work hrst with the value 3. As was already said, 
we have as an initial input a direct list which 3 is contained. For direct lists, 
we use an equation showed on 14: 

3.3! ^18 (16) 

After that, we observer the next element of the offset vector. As the previous 
element is in a direct list and it is odd, the element 0 will be in a inverted list. 
Using 15, we have: 

s^s + (2-0).2!^18 + 4^22 (17) 

Finally, observing the last element of the offset vector, we have the previous 
element which is in a inverted list and it is even. Looking at Table 3, we 
evidence that the element 1 is in a inverted list. Thus, we have: 

s ^ s +(1 - l).l! ^ 22 + 0 ^ 22. (18) 

After we have calculated the serial in a incremental mode, we add 1 to the 
serial, because the equations 14 and 15 depend on s — 1. Thus, the serial value 
characterizes it between (l,n!). Where n is the size of the permutation value 
given as an input: 

s ^ s + 1 ^ 22 + 1 ^ 23. (19) 

The serial 23 is the corresponding serial to the offset vector (1,0,3) which 
corresponds to the permutation vector (3,2,4,1). On the algorithm’s imple¬ 
mentation showed here, the increment of one unit to the serial was made in 
the beginning, before get into the loop. A boolean variable was specihed to 
determine if the list which belongs the elements is direct or inverted. The 
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conditions that determine if the list is direct or inverted is optimized from 
four to two conditions. The algorithm returns a non-negative integer value, 
corresponding to the serial number required. 

Serial Algorithm by Offset (SAO) 

Routine Specihcations: 

• s: Serial value; 

• n: Size of the offset vector; 

• i: Index of the offset vector; 

• d: Offset vector, alternating its elements on 0..n — 1; 

• direct: Boolean variable. It determines if the list is direct or not. 

s ^ [ dn-i.{n - 1)! ] -Fl 
direct <— true 

for i n — 1 to 1 do 

if {di = even and direct = true) or {di = odd and direct = false) then 
direct <— true 
s ^ s + di-ifi — 1)! 

else 

if {di = odd and direct = true) or {di = even and direct = false) then 
direct <— false 
s ^ s + {i - di-i).{i - 1)! 

end for 
return s. 


4 Serial Composition Method 

4-1 Basic Concepts 


On the combinatorial family, the composition of an integer n in k parts is 
dehned by: 


n = ri + r 2 + . .. + rk ri> 0 i = l..k 


( 20 ) 


Where the order of the elements is important on the compositions generations. 
The Next Composition Algorithm [1] does the task of making composition, 
obtaining the next composition starting from the one before, interactively 
working to reach until the last composition of the list. The proposal of the 
described method on this section is to obtain the vector from the list position. 
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with no need of processing the compositions one by one until getting the 
expecting vector. 


4-2 Construction of the SCM 


The SCM was built from the repetition patterns which are present on the 
composition vector throughout its list. It is possible to see that this repetition 
is done in specihc positions on the compositions lists: they can be obtained 
through calculus. This feature of the composition vector makes possible the 
cast among the specihc positions on the list. This way, it decreases the number 
of needing interactions to hnd the vector of a specihc position. The calculus 
of the specihc positions is dehned below. 


4-3 Used definitions on the built of SCM 


It is known that the total number of compositions of n in fc parts is dehned 

by 11]: 


J(n, k) = 


n + k — 1 
n 


( 21 ) 


It is also know that: 

k I 
\ 

\ 


n — % 


i=0 \ K — % 


n + 1 
k 


( 22 ) 


Associating 21 and 22 we have: 


^ln + k-2-i 
i=o \ n — i 


L 

^ n + k — l\ 


V ^ / 


(23) 


For a c composition vector of a /c length from the L composition list, the 
modihcation of kth element in L is determined by the index of the partial sum 
dehned on the equation 23. This way, it is possible to know when a component 
of c stops to repeat its current element to get modihcations. Taking as a basis 
equation 23, we can see this index modihcation with the dehnition of the M 
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matrix, exposed in 24. 


M(n, k) = 


^ n + k — 2 


n + fc — 3 
n 

n 

n 


n + k — ?) 

n — 1 

n + k — A 

n — 1 


n — 1 
n — 1 


V 0 ) 
^k-3\ 

[ 0 J 


(24) 


• Each line of M that corresponds to the elements of the sum on equation 23; 

• From one line to another, the superior element of the binomial is decreased 
in one unity; 

• M has k — 1 lines and n + 1 columns. 


Now let’s see bellow how M is used to represent the SCM execution. 


4-4 Description of the M Matrix Raster 


Initially, we have as initial information the serial number of the composition 
on the list, to obtain the composition vector. Let’s suppose s is the serial 
related to the input variable of the method to a composition vector with a k 
length. The elements of the first line of M are added until this sum is over 
s —1. When it’s found a value in the sum which is over s —1, the corresponding 
value of the column which the added element immediately before the actual 
one is found is given to the last element of the composition vector found. After 
finding the last element (the attributions of each element of the composition 
vector are done from the end to the beginning until the second element of the 
composition vector), the next element, that means the element before the last 
one of the composition vector is found in M going down in diagonal to the 
next line of M, re-starting the counting of indices starting from the actual 
column. This is done in a successive way, until the sum of the elements which 
were visited in M is s — 1. Let’s see an example of this strategy in a matrix 
M(7, 5) with serial s — 1 = 282, seen in 25. The result of the raster is the 
composition vector c(283) = {1,0, 2,1,3}. It is known that the first element 
of the vector does not need to be calculated using M, once it can be obtained 
from the complement of the sum of the elements already found through the 
raster done in M. There is in 25 on the first row with the partial sum until the 
third element (we allocate the number 3 on the last position of the list), the 
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second row with only one element of the partial sum (we allocate the number 
1 on the position before the last one of the list), and so on. The sum of the 
underlined numbers in 25 converged to s — 1. 


M(7,5) 


^120 mm 35 20 10 4 
36 28 21 15 10 6 3 1 
8 7654321 

^ 1 1111111 ^ 


(25) 


4-.5 Description of the SCM 


Related to implementation, the algorithm proposed here to the SCM is ab¬ 
stract to the structure of M matrix, dehned above. The loop is traced as 
dehned, applying on the hrst element of the composition vector the corre¬ 
sponding value of the complement of the sum of the elements yet to come. On 
the subroutine element there is the procedure to obtain each element of the 
composition vector. 

Serial Composition Method (SCM) 

Algorithm specihcations: 

• n: Number of the composition (Composition of n elements in k parts); 

• /c: Parts of the composition; 

• s: Serial of the composition; 

• a: Auxiliary variable which makes the increment to the convergence of the 
serial number; 

• x, y. Auxiliary variables, for the dehnition of new binomial indices; 

• Cij\ Combination of i elements j by j; 

• z: Value of the complement, used to apply the value to the hrst element of 
the composition vector. 

Routine: 

2 ; a ^ 0 
X ^ n + k — 2 
y 

for i ^ 0 to fc — 2 do 

Ck-i-i <— element 
Z ^ Z + Ck-l-i 

end for 
Co ^ n — z 

return c. 
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Subroutine element 
for j <— 0 to u — 1 do 

if a + Cx-i-j,y-j < s — 1 then 
d ^ a Cx—i—j^y—j 

else 

X <— X — j 

y ^y-j 

return j. 
end for 
return n. 


5 The inverse process of SCM (Unserial Method) 

Related to the inverse process of the SCM, the actual view is the composition 
serial number obtaining, having as input data the composition vector of a 
number n in k parts. In this case, each component of the composition vector 
is seen as part of superior intervals of a nested loops, which makes the raster 
on the composition vector. The inverse process has concise data input related 
to the data input of the SCM, once the inverse process will not work under 
the value convergence of a serial number method: everything is already well 
dehned on the composition vector components, known that is only necessary 
to process the loops with related interactions to each component of the com¬ 
position vector give as input. As to the SCM, the proposed algorithm for the 
inverse process to the SCM works only with indices referred to the dehned 
matrix on SCM. The use of the indices abstracting from the matrix structure 
decreases the memory use and processes only what is needed for the calculus 
of the trace to be run on M matrix. 

Serial Composition Method (Inverse process) 

Algorithm specihcations: 

• n: Number of the composition (Composition of n elements in k parts); 

• k: Parts of the composition; 

• s: Serial of the composition; 

• X, y: Auxiliary variables, for the dehnition of new binomial indices; 

• Cij\ Combination of i elements j by j. 

Routine: 

X u -|- /c — 2 
y 
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for i /c — 1 to 1 do 
for j ^ 0 to Cj — 1 do 

s < s + Cx-j-[[k-l)-i],y-j 

end for 

X j rp ft . 

' Jo L^2 

y ^y-Ci 

end for 
return s. 


6 Serial Partition of an n-Set Method 


6.1 Basic Concepts 


For the sets partitions, we consider a family of snbsets Ti, T 2 ,..., contained 
in a set S = {1, 2,..., n} that satisfy the conditions [1]: 


Ti n Tj = 0 (i # j) 

Ur. = s 

i=l 


(26) 

(27) 

(28) 


We don’t consider the order of the elements contained on each one of S snbsets. 
The algorithm Next Partition of an n-Set [1] hnds the next partition set from 
the cnrrent one, working only with local information. The section’s proposal 
is to obtain a partition of S related to its position on the list of partition, 
withont considering the information abont the partition vector. 


6.2 SPSM’s Construction 


The SPSM was bnilt from the analysis of the data ontpnt of the Next Partition 
of an n-Set algorithm. Each index that denotes the owned elements to a snbset 
has a pattern model and the search for these indices is made by a combinatorial 
strnctnre that take as basis a tree that represents all the partitions of a set 
with n elements, called Bell Tree. The nnmber of nodes of this strnctnre grows 
qnickly when n increases. So, the solntion for snch problem is to dehne a new 
strnctnre that does a mapping of the tree specihc positions. As a reference 
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we took a matrix structure to keep those positions, called D Matrix. The Bell 
Tree and D Matrix are dehned below. 


6.3 Bell Trees 


We can dehne a Bell Number as a number of partitions possibilities of a set 
with n elements. Such number is defined by [2]: 


B 


n 



(29) 


Where {2} is the Stirling number of second kind, which is defined by [2]: 



(30) 


With the definitions 29 and 30 above, this subsection proposes to show a 
combinatorial structure which is called here Bell Tree and the use of this 
structure to solve the Partition of an n-Set problem. It was used as basis the 
definition of a partition tree, shown in the algorithm Next Partition of an 
n-Set [1], associating the tree nodes with the equations which its terms are 
defined on the tree construction properties. 


6.4 Properties of Bell Trees 


Let n be the number of elements of a set that will be partitioned and B 
the Bell Numbers that denotes a set of terms that can form a mathematical 
expression in a node of the tree, the structure here dehned adopt the partition 
tree structure with n-Sets and added to it others properties. Follows below 
some properties of the structure: 

• The Tree has n levels; 

• The Tree nodes are non-negative integers numbers which is formed by ex¬ 
pressions that evolves exclusively Bell Numbers and multiplicative integer 
constants; 

• Let S' be a subset contained in B which denotes the number of distinct terms 
of an expression with a node N, we determine the number of descendants 
produced by N from the number of elements of S increased one unit; 
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Table 4 

Nodes values of the Tree illustrated on Figure 1 


Nodes 

Expressions 

1 


2 

Ba 

3 

Bq — Ba 

4 

Bs 

5 

Ba — B^ 

6 

Ba — B^ 

7 

Ba — B^ 

8 

(7^5 — Ba) — 2{Ba — B^) 


• The number R of the nodes on a Bell Tree corresponds to: 

n 

R = T.B, ( 31 ) 

i=i 

• Let k be the number of descendants of the actual node iV, the value for all 
its /c — 1 are defined by: 

Let be a mathematical expression with a w index and 2 ; terms from 
the ascendant node, the descendant node corresponds to To illus¬ 

trate such fact, if we have a ascendant node represented by the expression: 

Ew,Z = {Byj — Byj-l) — — B yj ^ 2 ) ( 32 ) 

Then, the descendants nodes will correspond to the expression: 

Ew-l,Z = {By,_l — By,_2) — 2{Byj_2 “ By,_^) (33) 

For the kth descendant of N, we have the expression: 

Ek,Z + l = Eyj^z ~ {k — l)Ey,_l^z^ (34) 

That it is equivalent to say that the kth descendant is equal to the as¬ 
cendant node value minus k — 1 times the value of the others descendant 
nodes. 

We can now make a tree that follows this formation law. For n = 5, we have a 
tree illustrated on Figure 1. The three first levels expressions of this tree are 
in Table 4, which the number of nodes are entered from top to bottom and 
from left to right. 


19 













Fig. 1. Bell Tree with level 5 


Ail 


6.5 Definition of D Matrix 

The D Matrix is a superior triangular matrix with n x n dimension. Its first 
column is made by the Bell Numbers 


(35) 



and the others columns are dehned by the equation below: 


(36) 



The matrix showed on 37 represents the D Matrix for n = 6. 


^203 151 77 26 6 1^ 

52 37 17 5 1 

15 10 4 1 

5 3 1 

2 1 





6.6 The D Matrix fulfilling 

Before process the SPSM, it is necessary hll in the components of D Matrix to 
speed up the attainment calculus of the search on the Bell Tree. The use of D 
Matrix abstracts the construction of the whole Bell Tree, processing only the 
necessary to the SPSM use. If the tree had been built completely, we would 
have the number of nodes R equivalent to the equation 31. Whereas on the D 
Matrix , we have 



( 38 ) 
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Algorithm for D Matrix fulfilling 

Algorithm Specifications: 

• D: Matrix that keep the specific positions of the Bell Tree; 

• i,j: Indices for the raster of D Matrix; 

• Bn'- Bell Number. 

Routine: 

for i 0 to n — 1 do 

Difl ^ Bji—i 

end for 

for i 1 to n — 1 do 
for j ^ 0 to n — i — 1 do 

Dj^i ^ Dj i—i i.D 

end for 
end for 
return D. 


6.1 SPSM’s Specifications 


After that the D Matrix has its values filled in, the SPSM can be invoked. We 
have an external loop that runs all the partition vector, attributing for each 
component by means of the sub-routine element that does a raster on the cur¬ 
rent tree level (the structure was abstracted from the method). On this raster, 
we have a condition that examines if the variable used for the indication of 
extrapolation surpassed the serial number token as an input. If it exceed this 
value, we go down a tree level; if not, we add the descendants’ current value to 
the control variable for later verification of extrapolation, token as an input. 
If the condition had not been satisfied for the whole loop, it means that the 
search on the tree’s level arrived to the last descendant of this one, indicating 
that the search will be expanded to the last descendant of the current tree’s 
level. 

Serial Partition of an n-Set Method 

Algorithm Specifications: 

• p: Partition vector; 

• D: Matrix that keep the specific positions of the tree; 

• i,j: Indices for the raster of D Matrix; 

• k'. Index for the raster on the p vector; 

• t'. Index that determines the element of each component of partition vector; 

• a: Number that does the convergence for the serial number given as input; 

• S'. Partition’s serial. 
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Routine: 


z ^ j ^ a <— 0 
for /c <— 0 to 77, — 1 do 
Pk ^ element 

end for 
return p. 

Sub-Routine element 
for t ^ 0 to j do 
if a -|- Di j > s then 
i ^ i + 1 

return t. 
else 

Qj *■— Cl Dij 

end for 
J ^ J + 1 

return j. 


6.8 Algorithm for stylized output of the partition vectors 

We know that the output gotten on both SPSM and Next Partition of an 
n-Set corresponds to the indices of each subset of the partition on the set. For 
the partition vector on position 26, we have the output (0,1,1, 0, 0) and the 
stylized output corresponds to (1,4, 5)(2, 3), indicating that 1,4 and 5 is on 
the hrst subset (indicated as 0), with 2 and 3 on the second subset (indicated 
as 1 on the partition vector). For show the stylized output, it was used a string 
vector to get the elements adequately on each subset which it belongs. Next, 
the algorithm does the output treatment and the all elements of each subset 
be ordered related to the last elements of partition as an input for the stylized 
output of the n-sets’ partitions method. 

Algorithm for stylized output data 

Algorithm Specihcations: 

• s: Strings vector that does the mapping of the elements referring to the 
partition vector’s indices; 

• p: Partition vector; 

• i: Indices for the raster of s vector; 

• r: Final output on string form; 

• n: Size of the partition vector; 

• +: Operator that denotes a concatenation between strings] 

• k: String that get the current element of the stylized output; 
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• length{w): Function that returns the size of the string w] 

• substring{w,ini, sup): Function that returns a w substring of the interval 
between ini and sup. 

Routine: 

for i ^ 0 to n — 1 do 

Spi ^ Sp^ + (* + !) + ”,” 

end for 

for i <— 0 to n — 1 do 

if length{si) > 1 then 

k ^ substring{si, 0, length{si) — 1) 
s, ^ ”(” + A: + ”)” 
if Si = ” 0” then 
r ^ r + Sj 
end for 
return r. 


7 The Inverted Process of SPSM (Unserial Method) 


For the inverted process of SPSM, we have a partition vector as input and the 
serial number as output. Such process is made taking each partition vector’s 
component as the number of loops that can does a sum of each specific position 
referring to the current node of the tree mapping through the D Matrix. 


7 .1 Inverted Process of SPSM Specifications 


Initially, we have the indices attribution referring to the D Matrix. The u in¬ 
dex has the attribution equals to 1 because it is not necessary do the mapping 
on the first component of the D Matrix, once it always begins with zero for 
being the main tree descendant. After that, we have an external loop that 
runs the whole partition vector given as input, with another loop that does 
the sum to the main positions’ serial mapped on D Matrix. After the sum at¬ 
tributions to the serial, we have the checking if the actual component should 
go down a line on D Matrix (equivalent to go down a level on the tree). If it 
should go down a level, the line is added by one unit. If not, the column on D 
is added by one unit, indicating the search to the descendants on the current 
node reached its last one and the search be doing on the other tree ramihcation. 

Inverted Process of SPSM 

Algorithm Specifications: 
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• p: Partition vector, given as input; 

• i,j: Indices referring to the raster on p vector; 

• D: Matrix of specific positions on the Bell Tree; 

• u,v: Indices referring to the mapping on D] 

• s: Partition’s serial, that is the result expected from de algorithm. 

Routine: 

i ^ j <— n <— 0 
s ^ M <— 1 

for i 1 to n — 1 do 
for jM— 0 to Pi — 1 do 

S ^ S ~\~ 

end for 

if j < Pi-i then 

M <— M + 1 

else 

n <— f + 1 
end for 
return s. 


8 Serial k-Subset of an n-Set Method 

8 .1 Basic Concepts 


We call (^) the number of possibilities to combine n things on k different parts. 
On analyzed literature [1], we have two ways to do this work on a sequential 
form. The hrst method builds the /c-subsets in a lexicographic order and the 
second method obtains the next subset from its predecessor, subtracting one 
element of the set and adding on other element of the subset. 

The proposal of this section is presents a method that, put in action iteratively, 
it creates a list in a lexicographic order given as input the position on the list 
of possible combinations. This strategy is more efficient when we desire to get 
a combination on a specihc position inside the whole list of combinations. 


8.2 Getting the combinations on lexicographic order 


The algorithm Next k-subset of an n-Set [1] is able to create on a simple way 
the combinations in lexicographic order in a non-recursive way. The recursive 
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model of this algorithm will be shown on this section optionally. The cur¬ 
rent combination of the recursive model is showed through the showOutput 
method, that is indicating here an output model of generic data for a com¬ 
bination that is going to have as output elements that can vary between 0 
and n — 1. The call to begin the routine must be done on a Combine{0) way, 
taking as stopped of the recurrence, the moment that the parameter did k 
recurrences. Next, we have a recursive model of the algorithm, without taking 
longer on its construction. 

Recursive /c-subset of an n-Set Algorithm 

Algorithm Specihcations: 

• k: Dimension of the subset; 

• s: Subset vector; 

• n: Set’s cardinality; 

• y. Vector that determines the changes of n set elements in subset k. 
Routine: 

Combine{i) 

for Si <— sum{i, i, 0) to n — {k — i) do 
if i 7 ^ fc — 1 then 
Combine{i + 1) 

else 

ShowOutput 
end for 
Vi+i ^ 0 

Pi ^ Pi+ 1 


Subroutine sum{w, j, z) 

for i ^ 0 to j do 

Z + Pi 

end for 
return z + w. 


8.3 SKSM’s Construction 


To build the SKSM, we observe the data output of the subsets’ list created 
by the Next k-subset of an n-Set algorithm, characterizing the repetition’s 
pattern in the elements of subset. In this case, the pattern can be delineated 
under a tree model. The formation law of this tree is showed below. 


25 


8-4 Definition of the Binomial Tree 


The structure that represents the repetition’s pattern is characterized as a 
tree formed exclusively by binomial coefficients. Let be the current node 
of the tree with w label, its descendants are dehned by 




( 39 ) 


/ n—k-\-l \ 

V fc-i ) 


w-\-n—k-\-l 


Each ascendant will have n + k — 1 descendants, and the repetition’s pattern of 
combinations are analyzed through the insertion of labels in each node of the 
tree. For each label of the current node, the descendants nodes enter its labels 
in relation to the ascendant node, indicating the value from each element of 
the subset founded on the tree. 


8.5 SKSM Specifications 

The SKSM abstracts the construction of the tree, doing the search calculating 
only the binomial coefficients and its relations with the nodes’ labels that were 
visited on the search. The method does an external loop attributing to each 
element of the subset the result of the subroutine element. The subroutine 
element does a search on the tree’s nodes. The abstraction of the tree is made 
through the indices changes of the binomial coefficient attributed on x and y. 
The method is described below. 

Serial /c-Subset of an n-Set Method 

Algorithm Specihcations: 

• p: fc-dimensional Subset; 

• n: Cardinality of the Set; 

• a: Auxiliary variable that is used to check the stopped of the method; 

• X, y: Indices of the binomial coefficients functioning as the formation law of 
the Binomial Tree; 

• r: Variable that controls the labels of Binomial Tree; 

• s: Serial of p subset; 
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Routine: 


a: <— n 
y ^k-l 
a <— r <— 0 

for i <— 0 to fc — 1 do 
Pi <— element 

end for 
return p. 


Subroutine element 
for j ^ 1 to a; — 1 / + 1 do 

< s then 


if a 

+ Cx-j,y 

a 

else 

<— a -|- U, 

X 

^ x-j 

y 

^y-1 

r 

^r+j 

return r. 


-j,v 


end for 
return r. 


9 The Inverted process of SKSM (Unserial Method) 


As many others serial methods from others combinatorial problems with same 
importance, the SKSM also has its inverted process. What is giving as input 
is the vector that represents the /c-dimensional subset with the cardinality of 
the set, called n, getting as output its position (i.e. the serial number) on the 
subsets list in a lexicographic order. 


9.1 SKSM Inverse Process Specifications 


On the inverted SKSM process, it is taken as the maximum limit of the in¬ 
ternal loop is the difference between the elements of the subset. This will 
determine how much the loop will repeat related to the elements of the sub¬ 
set. The indices of the binomial coefficient do an appropriate control for the 
adding among the nodes be made correctly. The inverted process also abstract 
the tree’s structure, getting the result by means of local information related 
to the subset given as input. 
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Serial fc-Subset of an n-Set Method (Inverted Process) 

Algorithm Specifications: 

• p\ fc-dimensional Subset; 

• n: Cardinality of the Set; 

• x, y. Indices of the Binomial Coefficients; 

• r: Variable that will control the internal loop of the method; 

• s: Serial of p subset; 

• c^,k- O- 


Routine: 

X ^ n 
y ^ k-l 
s <— 1 
r <— 0 

for i 0 to fc — 1 do 

for j ^ 1 to pj — r — 1 do 

S ^ S “ 1 “ Cx—j^y 

end for 

X ^ X — {pi — r) 

y 

r ^ Pi 

end for 
return s. 


10 Serial Subset Method 

10.1 Basic Concepts 


We have 2” possible configurations to get the subsets of a {1, 2,, n} set. The 
disposition of the subset’s elements can be represented by a flag which acti¬ 
vates its insertion into the subset. For example, the configuration {1, 0,1,1, 0} 
represents the subset {1, 3,4}. The Next Subset of an n-Set [1] algorithm does 
this job in a sequential form, returning the next subset from the actual one. 

This section proposes to present a method that returns a subset of an n-Set 
from its serial number. The inverted process for this method is also showed, 
getting the serial number from the subset given as input. 


10.2 SSM’s Construction 


The SSM was built from the analysis of repetition’s pattern that the subsets 
presents in the complete list of subsets. Similarly to the Serial Permutation 
Method related to repetition’s pattern of the offset vector, the SSM has a 
regular repetition pattern, being able to get each subset’s component through 
a closed equation. 


Let p be a subset of a n-Set, each subset’s component with index /c = 0,1,..., n— 
1 related to a serial number s = 1, 2,..., 2"' is defined by 


Pk 


(s — 1 + 2^) mod 2^+^ 

2^+i 


( 40 ) 


The equation 40 doesn’t present restrictions in the pattern’s representation. 
Each component pk € {0,1}, and 2^ in (s — 1 + 2^) is related with a cyclic 
structure which the numbers are presented at the list. This peculiarity is easily 
identified on output data of Next Subset of an n-Set [1] algorithm and the 
method here proposed. 


10.3 SSM Specifications 


About the SSM’s specifications, we have a loop that do an association between 
the equation 40 and the structure that represents the subset. The method re¬ 
turns a subset from the serial number given as input, where the number 1 
on the output data indicates which subset’s component will be activated, as 
explained in the basic concepts on this section. 

Serial Subset Method 

Algorithm Specihcations: 

• p: Subset of n-Set; 

• n: Cardinality of the Set; 

• s: Serial of p subset; 

Routine: 

for i ^ 0 to n — 1 do 

Pi ^ L((s -1 + 2*) mod 2*+2) / 2*+^J 

end for 
return p. 
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11 The Inverted Method of SSM (Unserial Method) 


In this section the inverse process of the SSM will be shown. The input data 
is the subset (the elements of the set that will be activated) and the output 
data is the corresponding position in the list of subsets. In the case at hands, 
we determine which it is the moment to walk through specihc positions of the 
identified repetition’s pattern in subsets. The inverted process for subsets also 
have similarities with the inverted process of Serial Permutation Method. The 
list of subsets, from right to left, presents a sequence of O’s and I’s, in this 
order. From there, the list for analysis of the subset’s next component can 
be inverted in the case of this subset’s element be equals to 1. In this case, 
we cannot advance in the list for convergence of the serial’s desired value. 
Otherwise, the serial value is modihed on iterative form among the specihc 
position of the list until reach the subset’s serial. 


11.1 SSM Unserial Method Specifications 


In the SSM’s specihcation, it was dehned a logical variable which dehnes if 
the list of current subset’s component is direct or inverted, dehning as true for 
direct lists and false for inverted lists, starting this variable as true because 
the components’ disposition has been analysed from right to left. The initial 
serial’s value is 1. After that, we have a decreased loop that does the verihca- 
tions in each component of the subset for determines if it advances or not for 
the convergence of the hnal result. It was observed that such condition in the 
way as was constructed can be represented by eXclusive OR logical connec¬ 
tive, usually denoted by xor operator. When the xor is satished, the logical 
variable is activated as true, indicating that the next component is included 
on a direct list. Otherwise, the serial advances other specihc position and the 
list is conhgured as inverted. 

Inverted Process of the SSM (Unserial Method) 

Algorithm Specihcations: 

• p: Subset of n-Set; 

• n: Cardinality of the Set; 

• s: Serial of p subset; 

• d: Logical variable, which indicates if the list of subsets is on direct or 
inverted order. 

Routine: 

s 1 
d <— true 
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for i ^ n — 1 to 0 do 

if {Pi = 1 xor d) = true then 
d ^ true 
else 

d ^ false 
s ^ s + 2* 

end for 
return s. 


References 

[1] NIJENHUIS, A., WILF, Herbert S. Combinatorial Algorithms for computers 
and calculators. Academic Press, INC, 1978. 

[2] GRAHAM, R. L., KNUTH, D. E., PATASHNIK, O. Concrete Mathematics: 
A Foundation for Computer Science, second edition, Addison-Weslay, Reading, 
Mass., 1994. 

[3] WILF, H. S. East Side, West Side, online version, 2000, 
http://www.cis.upenn.edu~wilf/. 


31 



